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Abstract 



Using the angular correlation between the vr^ emitted in a D*^ — > 
D^ir^ decay and the e"*" emitted in the subsequent — > Xe^iy de- 
cay, we have measured the branching fraction for the inclusive semi- 
electronic decay of the meson to be: 



B{D^ Xe+iy) = [6.64 ± 0.18{stat.) ± 0.29{syst.) 



70. 



The result is based on 1.7 fb~^ of e^e~ collisions recorded by the 
CLEO II detector located at the Cornell Electron Storage Ring (CESR). 
Combining the analysis presented in this paper with previous CLEO 
results we find, 

^ = 1.684 ± 0.056(stat.) ± 0.093(sysi.) 

B\D" — ^ K TT'^) 

and 

B{D K-e+u 
B{D Xe+v) 



0.581 ± 0.023(stat.) ± 0.028(syst.). 



The difference between the inclusive rate and the sum of the measured 
exclusive branching fractions (measured at CLEO and other experi- 
ments) is (3.3 lb 7.2)% of the inclusive rate. 



1 Introduction 

In this paper, we present a new measurement of the inclusive semi-electronic 
branching fraction of the meson. The comparison of the measured inclu- 
sive semi-leptonic branching fraction with the sum of the observed exclusive 
semi-leptonic branching fraction provides a measure of missing or unobserved 
modes. Recent experimental progress on exclusive measurements has yielded 
precise measurements of the dominant Cabibbo favored modes, observation 
and measurement of the Cabibbo suppressed branching fractions and strin- 
gent upper limits on suppressed Cabibbo favored branching fractions, but has 
not yielded an improvement in the measured inclusive semi-leptonic branch- 
ing fraction. The inclusive branching fraction measurement presented here 
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and previous measurements of exclusive branching fractions allows for more 
accurate comparison than previously performed. For a complete review of 
experimental and theoretical developments we refer the reader to recent re- 
views ii. 

In addition, we combine the inclusive result presented here with previous 
CLEO results on B{D^ R-ir) and B{D^ K-e+v)/B{D^ K-ti+) 
[|, I to obtain the ratio, B{D^ K-e+v)/B{D^ Xe+z/). As a check 
of the method the observed inclusive electron momentum spectrum is also 
extracted from the data and compared with a Monte Carlo simulation. 

2 Analysis Technique and Event Selection 

The technique to measure the absolute inclusive semi-electronic branching 
fraction of mesons is similar to the previous CLEO absolute branching 
fraction measurement of K^tt^ P]. Common to both analyses is the 

method used to determine the number of D*^ D^n^ decays in the data 
with minimal systematic bias. It is based on the unique two body kinematics 
of the D*~^ —>■ D^TT'^ decay and the topology of e^e~ — cc reactions at a 
center of mass energy of 10.5 GeV. Briefly, the idea is that the thrust axis 
(defined to be that axis along which the projected momentum is a maximum) 
for the event approximates the D*~^ direction in the lab. The limited amount 
of available phase space in the D*^ D^n^ decay, results in a small angle 
between the thrust axis and the charged pion. We denote this angle between 
the thrust axis and the charged pion as a. Also, the magnitude of the 
pion momentum is correlated to the parent D*~^ momentum. Pions with 
momentum greater than 225 MeV/c are kinematically forbidden to come 
from the T(4^) BB, B D*+X, D*+ D%+ decay chain. This 
selection assures that the D*'^ is from e+e^ cc production and the event 
has a well defined thrust axis. The top plot in Figure |I| shows the sin^ a 
distribution for all pions with momentum between 225 and 425 MeV/c in 
the data. The peaking at low sin^ a is evidence for D*^ D^-n^ decays. 
The total number, N{D*^ D^n^), of decays in the sample is 165658 ± 
1149(stat.) ±2485(syst.). This total is identical to that presented in Ref. 
as the same data and selection criteria are used in both analyses. 

The total number of semi-electronic decays, N{D*'^ —>■ D^n^, 
Xe'^u), is determined by identifying an within a cone around the vr"'" di- 
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rection and plotting the sin^ a distribution for those tt^ with an associated 
. This is achieved by studying the sign correlated vre combinations in the 
data. "Right sign" combinations, vr+e^, provide the signal distribution and 
"wrong sign" combinations, vr+e", are studied to aid background determi- 
nations. Once the number of Xe^v decays has been determined, the 
branching fraction is then. 



where e{D^ — > Xe^v) is the efficiency for detecting the electron. 

A detailed description of the CLEO II detector can be found in Ref. p. 
Electrons and positrons are identified principally from the ratio of the en- 
ergy measured by the Csl calorimeter and the momentum measured by the 
drift chamber (E/p). Additional information on energy loss in the drift cham- 
ber and shower shape in the calorimeter is also used to maximize the identifi- 
cation efficiency and minimize the mis-identification of hadronic tracks. The 
electrons are required to have momentum greater than 0.7 GeV/c and a polar 
angle with respect to the beam axis {9) between 45° and 135°, to insure a well 
determined efficiency and minimal uncertainty due to mis-identified hadronic 
tracks. Furthermore it is important to reduce the number of electrons from 
Xvr^; 7r° e"'"e~7, where the e"'"e~7 final state is due to either a 
Dalitz decay of the vr*^ or a 7 conversion in the detector material. This is 
accomplished by requiring that the identified electron, when combined with 
each opposite sign track in the event, does not yield an electron-positron 
mass below 0.050 GeV/c^. Every opposite sign track is used to form these 
pairs, whether or not it is identified as an electron. 

In order to correlate the vr"*" with an a fiducial angle cut is applied in 
the lab frame. We require that cos(9e-7r) > 0.8, where Qe--K is the angle 
between the vr emitted in the initial D*^ D^tt^ decay and the electron 
from the subsequent —>■ Xe^v decay. The bottom histogram in Figure |l| 
shows the sin^a distributions for vr's after requiring an electron within this 
angular region; the solid squares are for ix^e^ combinations (right sign) and 
the open squares are for Ti^e~ combinations(wrong sign). 



B{D' 



,0 



Xe+i/) 



N{D*+ D\+, Xe+v) 




N{D*+ DH+) X e(L)0 ^ Xe+u) 
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Figure 1: The inclusive sin^a distribution for candidate pions (open circles) 
and the derived non-D*'^ background (solid line) in the top plot. Requiring 
an electron near the pion with the same (opposite) sign results in the solid 
(open) squares in the bottom plot. 
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3 Extraction of yields 



As previously stated, the yield of D*^ D^n^ decays is identical to that 
presented in Ref. P]. In this section we detail the determination of the 
number of —> Xe^v decays associated with the initial D*^ D^tt~^ 
decay. 

The sin^ a distribution for vr+e^ (right sign) combinations contains three 
distinct components: signal and two types of background. One background 
has a sin^ a distribution that is identical to the signal as it originates from the 
decay D*~^ — > D^tt"*", X/^+e+, where fn+e+ denotes a e"*" from either a 

hadronic track mis-identified as a electron or an electron from a vr" — > e~^e~'y 
final state. The other background is due to random soft pions (225 to 425 
MeV/c in momentum), in coincidence with an electron, and is not as sharply 
peaked near sin^ « = as the signal distribution. 

The sin^ a distribution for 7i~^e~ (wrong sign) combinations is devoid of 
signal but contains the same two sources of background as the right sign 
distribution 0. The shapes for these backgrounds in the right sign and 
wrong sign distributions are identical, although the normalizations differ. 
This difference in normalization is the result of the hadronic track mixture 
{ir/K ratio) combined with the hadronic track mis-identification rates. For 
the non-Z)*"*" pion background, the normalization is different due to charge 
conservation in the event. 

To use as much information as possible, the right sign and wrong sign 
distributions are fit simultaneously to the following functional forms: 

grs{P7v,sm'^ a) = NDO^xe+APn)gDO^Xe+uism'^ a,p^) + 

Brsip.)P2ism^ a) (2) 

g^s (Ptt, sin^ a) = Nr,o_x {P7t)9do^x (sin^ a,p^) + 

-B^s(p7r)-P2(sin^a). (3) 

The expected sin^ a distributions, (7£)0_,xe+v(sin^ a,Pn) and (7£)0^x/e(sin^ a, 
for the right sign and wrong sign distributions are obtained from a Monte 
Carlo simulation. The Monte Carlo simulation correctly reproduces the mea- 
sured D*~^ production momentum distribution, and simulates — > Xe^u 
decays via the "cocktail" of exclusive modes presented in Appendix A. The 
second order polynomial, P2, is constrained to have the same shape for 
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p{7r) 




Yields 


(MeV/c) 


Right Sign 


Wrong Sign 


225 - 


250 


1232 ± 53 


32 ±31 


250 - 


275 


1071 ±49 


74 ±29 


275 - 


300 


935 ± 44 


45 ±25 


300 - 


325 


689 ± 38 


39 ±22 


325 - 


350 


414 ± 32 


-29± 18 


350 - 


375 


259 ± 25 


36± 17 


375 - 


400 


166 ±20 


-4 ±12 


400 - 


425 


79 ±15 


± 11 


Total 




4845 ± 104 


193 ± 62 



Table 1: The total yield of right sign and wrong sign events as a function of 
pion momentum. Backgrounds have not yet been subtracted. 



both the wrong and right sign distributions. The yield of D° Xe'^u 
decays, NDO_xe+u{Pn), the yield of mis-identified hadrons and electrons from 
7r° — i> e^e~7 in the wrong sign distribution, N]:,o^xfXP-^)i ^^e normalizations 
and shape, Brs{pn), B^sip-w) and P2, of the background polynomial are de- 
termined from the fits to the sin^ a distribution of the right sign and wrong 
sign samples in bins of pion momentum, p^. 

The sin^ a distributions for the data, with the resulting fits overlaid, are 
shown in Figures ^ and |^ and the right sign and wrong sign yields are pre- 
sented in Table |l|. The right sign yields reported in this table have a contri- 
bution due to X/jr+e+ backgrounds. 

4 Determination of the background contribu- 
tion to the signal 

In this section the magnitude of the /7r+e+ background to the right sign 
signal yield is determined. The two contributions to this background are 
the following decay chains: D*^ D^tt~^, —>■ Xh^, where the h'^ is a 
hadronic track mis-identified as an electron and D*'^ D^n'^, Xtt", 
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Figure 2: The sin^ a distribution for pions with momentum between 225 and 
325 MeV/c with an identified electron within cosOtt-e > 0.8. Events with 
the electron and pion having the same sign (right sign) are plotted on the 
left side, the opposite sign events (wrong sign) are plotted on the right side. 
The points represent the data and the histogram is the result of the fit. The 
dashed line represents the random pion-electron background and is modeled 
by a second order polynomial. 
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Figure 3: The sin^ a distribution for pions with momentum between 325 and 
425 MeV/c with an identified electron within cos07r_e > 0.8. Events with 
the electron and pion having the same sign (right sign) are plotted between 
on the left side, the opposite sign events (wrong sign) are plotted between on 
the right side. The points represent the data and the histogram is the result 
of the fit. The dashed line represents the random pion-electron background 
and is modeled by a second order polynomial. 
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TT 



.0 



e+e 7. We denote this sum for the right sign background as 



rs 



iV(X7r°)^+e+(X7rO) + NiXh+)U,^{Xh+), 



(4) 



where N{X7c^) [N{Xh+)] is the number of inclusive D° ^ Xn^ [D^ Xh+] 
decays in the data and f.n-+e+ (^7r°) [f-K+e+ i^h'^)] is the efficiency for detecting 
this background as signaL We can define the same sum for the wrong sign 
yield as 



The only difference between the wrong sign yield and the right sign back- 
ground is due to the fact that the positive tracks from decays are much less 
likely to be kaons than negative tracks from decays. Using fj^+e- {Xtt^) = 
f^+e+{X7[^), we find 



Nrs = N^s-N{Xh-)U+,-{Xh-) + N{Xh^)U+,+ iXh^). (6) 
If N{Xh+)f^+e+{Xh+) = N{Xh~)f^+e-{Xh-), then the wrong sign yield 



would be equal to the background contribution to the right sign yield. How- 
ever, the TT"*" : ratio of tracks originating from mesons is quite 
different from the vr" : K~ ratio. Using world averages of the measured 
branching fractions, the 7V~^:K~^ ratio is 96:4 while the n~:K~ ratio is 
42:58, for pions and kaons from mesons which pass the same geometry 
and momentum criteria as for electrons. This difference coupled with differ- 
ent mis-identification rates for vr's and K's leads to a small correction to the 
wrong sign yield. 

The probability for a 7r+ track to be mis-identified as a e+ is determined 
by studying a large sample tt+tt^ decays in the data. This sample is 

large enough to determine the mis-identification probability for charged pions 
as a function of their momentum. This probability is measured to be (0.056 ± 
0.015)% for pions with momentum between 0.7 and 0.9 GeV/c. It rises as 
a function of pion momentum, such that for pions with momentum between 
1.9 and 2.5 GeV/c, it is measured to be (0.250 ± 0.059)%. Convoluting 
the momentum dependent mis-identification probability with a Monte Carlo 
simulation of the 7r+ momentum distribution from and decays, we 
find the mis-identification probability integrated over pion momentum to be 
(0.102 ±0.016)% for the right sign pions and (0.093 ± 0.011)% for the wrong 



ws 



N{X7r')U+,-{X7r') + NiXh-)U^,-iXh-). 



(5) 
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sign pions. These numbers differ due to the different momentum spectrum for 
right sign and wrong sign pions. The error is due to the statistical uncertainty 
in the mis-identification probabihty per track as a function of momentum. 

For charged i^'s the data do not provide a statistically rich and clean 
sample as for pions. The cleanest sample of charged kaons comes from recon- 
structed ir-7r+(7r°) decays. With 19742 ±221 reconstructed D^'s with 
a K~ that passed the momentum cuts, 4.5 ±5.5 were consistent with the K~ 
being identified as an electron. This yields a central value of (0.02 ± 0.03)% 
for the mis-identification probability due to kaons. As no momentum depen- 
dence measurement is possible we use (0.02 ±0.03)% as the mis- identification 
probability for charged kaons over the whole momentum range of interest. 

Multiplying these mis-identification probabilities by the vr : K fractions, 
we find the total mis-identification probability of fj^+e+iXh^) = (0.099 ± 
0.016)% for the right sign hadronic tracks and /^+e-(X/i^) = (0.051 ± 
0.016)% for the wrong sign hadronic tracks. This is almost a factor of two 
difference between the two mis-identification rates. These mis-identification 
probabilities represent the rate per hadronic track from D*^ D^tt^ de- 
cays where the vr"*" had momentum between 225 and 425 MeV/c. Since 
the extraction of yields is done in eight 25 MeV/c momentum bins, these 
mis- identification probabilities are determined for each of the eight bins indi- 
vidually. Small variations arise due to different momentum spectra and 
small changes in the vr : K ratio. 

To turn these mis-identification probabilities into the actual yield of mis- 
identified tracks, the inclusive right sign and wrong sign rate [N{Xh'^) and 
N{Xh~)] is determined from the data. The number of wrong sign and right 
sign hadronic tracks associated with D*^ —>■ D^7r~^ decays is determined 
by using the same code and technique as for identified electrons, with the 
requirement that the hadronic track not be identified as an electron. The 
resulting estimated mis-identified charged track contribution to the right and 
wrong sign yields is given in Table ^ as well as the final estimated total 
background to the right sign yield. 

5 Efficiency 

The efficiency for detecting the e"*" determined by the Monte Carlo simula- 
tion depends on the cocktail of exclusive modes used to generate the inclusive 
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p{7l) 

(MeV/c) 






Fh+ 


-'■"rs 


225 - 


250 


32 ±31 


13 ±3 


24 ±3 


43 ±31 


250 - 


275 


74 ±29 


13 ±3 


22 ±3 


83 ±29 


275 - 


300 


45 ±25 


9±2 


17 ±2 


53 ±25 


300 - 


325 


39 ±22 


7±2 


13 ±2 


45 ±22 


325 - 


350 


-29± 18 


5± 1 


9± 1 


-25 ± 18 


350 - 


375 


36 ±17 


4± 1 


6± 1 


38± 17 


375 - 


400 


-4± 12 


2± 1 


4± 1 


-2 ± 12 


400 - 


425 


0± 11 


1 ± 1 


2± 1 


1 ± 11 


Total 




193 ± 62 


54 ±5 


97 ±6 


236 ± 64 



Table 2: Summary of the expected background contribution as a function 
of pion momentum to the right sign yield, where A''^^ = N^s — F^- + Fh+ , 
Fh- = N{Xh-)U+,-{Xh~) and Fj,+ = N{Xh+)U+,+ {Xh+). 

semi-electronic decays. The ratios of exclusive rates presented in Appendix 
A are used to calculate the ratios Xm = B{D^ me"*"//)/ J2n B{D^ ne'^u), 
where m,n = K' , K*' , K{ (1270), K*' {U30),7r' , and p" mesons. The effi- 
ciency for each of these modes is obtained from a Monte Carlo simulation of 
each individual mode. The inclusive efficiency is obtained from 

e(Xe+z/) = X e(L'° ^me+z/). (7) 

m 

The extraction of yields is done in eight pion momentum bins from 225 to 
425 MeV/c, as in the K~7r^ analysis. Table ^ contains the efficiency 

in each of the eight pion momentum bins. The efficiencies for the individual 
exclusive channels are in Table |l^ (Appendix A). The total systematic error 
due to uncertainties in the cocktail is determined by varying the ratios in 
Table by one standard deviation, individually and collectively. The largest 
variation in the overall efficiency is seen when Xk and X^r are both raised 
or both lowered and the other modes are changed in the opposite direction. 
This causes a ±2% change in the efficiency and is the estimated systematic 
error due to the uncertainties in the cocktail of exclusive modes. 

In addition to changing the cocktail the effect of the assumed depen- 
dence of the form factors is studied by changing the ISGW slope (k) 
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The value used to generate the decays is k = 0.57 ±0.07, measured in a large 
sample of K~e^v decays by CLEO 0|. Variations of one sigma on k 

resulted in a ±0.6% variation in efficiency. The longitudinal and transverse 
contributions from K*~e^v decays were varied by one sigma of their 

measured value and the total efficiency changed by less than ±0.08% 0. 

6 Results 

6.1 B{D^ Xe+v) 

The relevant measurements for determining B{D^ Xe^v) are given in 
Table |^. The first column gives the inclusive D*^ D^ti^ yields from 
Ref. 0, the second gives the background subtracted yield of Xe^u 
decays, followed by a column of efficiencies. The last column is the branching 
fraction for Xe^v for the eight momentum bins. As a check that 

the eight measurements are self consistent, the was calculated under the 
assumption that all eight branching fraction measurements come from the 
weighted average. The result is a of 9.4 for 7 degrees of freedom. 

Sources of systematic effects and their estimated magnitude are listed 
in Table ^. The dominant systematic uncertainty is the evaluation of the 
electron identification efficiency. The electron identification algorithm was 
developed using clean radiative Bhabha events in the data sample. Its per- 
formance on continuum events is studied using vr" ^e^e~ where the e^e~ 
pair could originate from either a Dalitz decay of the 7r° or a 7 conversion 
in material. This study resulted in a conservative estimate of the electron 
identification systematic uncertainty of ±3%. 

The inclusive semi-electronic branching fraction is measured to be 

Xe+z/) = [6.64 ± 0.18 ± 0.29]% (8) 

where the first error is statistical and the second error is the estimated system- 
atic effect. Sources of model dependence have been minimized by relying on 
the experimental measurements of the exclusive rates of the observed modes 
and experimental measurements of the dV/dq^ spectrum in K'^e^u 
decays. Models have been used only for the dV/dq^ spectrum of the other 
exclusive modes. The previous value of [7.01 ± 0.62]% agrees with this re- 
suh §. 
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p(7r) 


N(D*+ D°7r+) 


N{D* 




e(Xe+z/) 


B(D^ Xe+z/) 


(MeV/c) 




D' 


° ^ Xe+z/) 


(%) 


(%) 


225 - 250 


44161 ±611 




1189 ±61 


37.9 


7.10 ±0.38 


250 - 275 


39114 ±562 




988 ± 57 


40.1 


6.30 ±0.38 


275 - 300 


29482 ± 475 




882 ± 51 


42.7 


7.01 ± 0.42 


300 - 325 


21120 ±396 




644 ± 44 


43.7 


6.97 ±0.49 


325 - 350 


14973 ± 334 




439 ± 37 


45.5 


6.42 ±0.56 


350 - 375 


9165 ±267 




221 ± 30 


48.0 


5.02 ±0.70 


375 - 400 


5492 ± 208 




168 ± 23 


49.5 


6.18 ±0.88 


400 - 425 


2151 ± 147 




78± 19 


50.7 


7.15 ±1.8 


Total 


165658 ± 1149 




4609 ± 121 




6.64 ±0.18 



Table 3: The yields of inclusive D*^ —>■ Z}°7r+ decays, Xe+z/ decays, 

the efficiency for detecting the Xe+i/ final state and the calculated branching 
fraction, as a function of the initial /?*+ pion momentum. The errors are sta- 
tistical only and include the statistical error on the background subtraction. 



Source Estimated Systematic Error 

(%) 

Electron Identification Efficiency ±3.0 

Xe+z/ Cocktail ±2.0 

N{D*+) ±1.5 

Track Reconstruction ±1.0 

Monte Carlo Statistics ±1.0 

Electron fake rate ±1.0 

Form Factor slope k ±0.6 

Total ±4.3 

Table 4: Estimate of the systematic uncertainty in the measurement of 
B(L>o ^ Xe+iy) 
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p(7r) 


N{D*+ 




e{KTi) 


N{D*+ D^7i+, 


e{Xe+u) 




(MeV/c) 






(%) 




(%) 


225 - 250 




1129 ±44 


64.6 


1189 ±61 


37.9 


1.80 ±0.12 


250 - 275 




945 ± 40 


64.3 


988 ± 57 


40.1 


1.68 ±0.12 


275 - 300 




741 ± 34 


64.4 


882 ±51 


42.7 


1.80 ±0.13 


300 - 325 




528 ± 30 


65.1 


644 ± 44 


43.7 


1.82 ±0.16 


325 - 350 




393 ± 25 


66.0 


439 ± 37 


45.5 


1.62 ±0.18 


350 - 375 




262 ± 19 


66.4 


221 ±30 


48.0 


1.17±0.18 


375 - 400 




153 ±15 


68.8 


168 ± 23 


49.5 


1.53 ±0.26 


400 - 425 




57 ±9 


63.1 


78± 19 


50.7 


1.70 ±0.50 


Total 




4208 ± 83 




4609 ± 121 




1.684 ±0.056 



Table 5: The yields of D ^ K^n^ decays, the efficiency for detecting the 
K~7T^ final state, the yields of D° — ^ Xe+z/ decays, the efficiency for detect- 
ing the Xe'^u final state and the calculated ratio of branching fractions, as 
a function of the initial D*~^ pion momentum. The errors on the data yields 
are statistical only. The ratio of branching fractions error is statistical only. 



6.2 B{D^ Xe+iy) /B{D^ K-7T+) 

In addition to measuring the absolute Xe^v branching fraction, it is 

straightforward to combine the yields presented here with those in Ref. 
to obtain a measurement of the ratio B{D^ —>■ Xe^v)/B{D^ —>■ K~7i^). 
This tabulation is done in Table This ratio is independent of systematics 
associated with the inclusive D*^ D^n^ yields. The contributions to the 
systematic error are given in Table ^. The result is 

^ Xe-^u)/B{D° K-'K+) = 1.684 ± 0.056 ± 0.093. (9) 

Again the first error is statistical and the second error is the estimated sys- 
tematic effect, where the use of a common dataset allowed cancelation of 
some systematic effects present in the individual results. 
This ratio allows for a check of the ratio 

Xk = B{D ^ K~e+u)/B{D Xe+u) (10) 
B{D K-e+v) B{D K-71+) 
B{D K-n+) ^ B{D Xe+u) ^ ^ 
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Source Estimated Systematic Error 

TO 

Electron Efficiency ±3.0 

Xe+u Cocktail ±2.0 

Track Reconstruction ±3.8 

Monte Carlo Statistics ±1.2 

Electron fake rate ±1.0 

Form Factor slope k ±0.6 

K~7T^ (mass fit and momentum cut) ±0.7 

Total ±5.5 

Table 6: Estimate of the systematic uncertainty in the measurement of 
^ Xe+u)/B{D^ K-71+). 



which is used in the — s- Xe^u cocktail. To obtain the most precise 
value possible, we take advantage of the fact that the CLEO results for 
B{D^ K~e^v)/B{D^ K~7c~^) were obtained with the same detector, 
allowing reduction in the systematic bias due to lepton identification (re- 
duced to ±1.7%) and the systematic bias due to tracking reconstruction 
(reduced to ±2%). There is also a large overlap of K~tt^ events 

which were used to calculate the two ratios which appear in Eq. ^ [|I^]. 
Using only CLEO results and taking these common systematic effects into 
account we obtain X^ = 0.581 ± 0.023 ± 0.028. Using all measurements 
of B{D^ — > K~e^v)/B{D^ K^tt^) and taking advantage of the common 



CLEO systematic errors results in a value of X^ = 0.552 ±0.035 [|1^]. These 



results agree well with the input value of Xk listed in Table 12 



6.3 Comparison of inclusive measurement to the sum 
of the exclusive rates 

The measurement of the inclusive semi-electronic branching fraction is often 
compared to the sum of the measured exclusive channels |T^. This compari- 



son provides a measure of the consistency of the experimental measurements. 

In terms of the branching fraction ratios, Rm = B{D^ — me^v)/B{D^ 
K~e^v), which are used in Appendix A for tabulating the —>■ Xe^u 



19 



cocktail listed in Table p!2l the ratio between the difference of the inclusive 



and the sum of the exclusive rates can be written as: 



B{D^ Xe+v) 



1 — Xi^(\ + Rk* + R-iT + Rpj- 

(12) 



Performing the comparison using only CLEO data {Xk = 0.581 ± 0.036, 
and 1 + Rk* + R-k = 1-724 ± 0.078) resuhs in a value of: 

^ xe+u) - (-0-2 ±7.7)%. (13) 

This CLEO result does not include a contribution from Rp as CLEO has not 
reported a value for this ratio. Inclusion of the small contribution for Rp will 
result in a central value further from zero, while still entirely consistent with 
zero given the experimental errors. Using the value of Xk = 0.552 ± 0.035 
obtained in the previous section and 1 + Rk* + Rn + Rp = 1.751 ± 0.067 (see 
Table |g) we find, 

^(^0 _ xe^u) - (3.3 ±7.2)%. (14) 

These results are consistent with the upper limits obtained by direct searches 



for the unobserved exclusive modes |T3|. 



6.4 The inclusive electron momentum spectrum 

The lepton spectrum from semi-leptonic charm decays has not been updated 



since the DELCO results |T^. Because the measurement presented here is not 
made in the rest frame of the we compare the observed lepton spectrum 
in the lab frame with that of the Monte Carlo simulation. To obtain the 
momentum spectrum for inclusive D*^ Xe^v decays, events were selected 
if they pass all the selection criteria previously described. An additional 
cut of sin^ a < 0.12 is applied. This cut retains 90% of the signal and is 
large enough that systematics associated with modeling the thrust axis are 
minimized. There is still background in this sample whose shape is provided 
by the wrong sign sin^ a distribution. The normalization of this background 



is obtained by normalizing the wrong sign sin^ a distribution to the right 
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sign sin^ a distribution for values of sin^ a > 0.2. As this result is focused 
on the distribution of the electron momentum not the normalization, a ±1% 
uncertainty in the level of background spread over the momentum range is 
negligible. In Figure ^ the background subtracted momentum spectrum for 
the electrons is shown along with the momentum spectrum obtained from 
the Monte Carlo simulation. The two distributions are normalized to the 
same number of events, resulting in a 75% confidence level. The comparison 
shows that the simulation is correctly producing D*^, mesons and the 
inclusive — > Xe'^u decays. Any deviations would indicate a problem in 
the simulation, either in the production or decay dynamics. We conclude 
that the Monte Carlo provides a good simulation of the data. 

7 Conclusions 

We have presented a new measurement of the inclusive branching fraction 
for — i> Xe'^u decays. The final result is, 

B{D^^Xe^v) = [6.6A±0.18{stat.) ±0.29{syst.)]%. (15) 

We find that the difference between this inclusive rate and the sum of the 
observed exclusive channels is (3.3 ± 7.2)% of the inclusive rate. This corre- 
sponds to an upper limit on the unobserved modes of 14% of the inclusive 
rate (at the 90% C.L.). The experimental upper hmits obtained using direct 
searches for specific unobserved exclusive semi-electronic modes are lower 
than the limit quoted here. However, the upper limit obtain in this paper is 
less sensitive to the assumption of what exclusive channels are unobserved. 
The two methods, direct searches and inclusive-exclusive rate comparison, 
both suggest that the remaining unobserved exclusive semi-leptonic modes 
occur at small rates. In addition the observed electron momentum spec- 
trum from inclusive Xe^v decays is seen to be well described by the 
exclusive semi-electronic cocktail. 
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Figure 4: The lab momentum spectrum of electrons from semi-electronic D° 
decays. The sohd squares represent the background subtracted data, and the 
histogram is the result of a Monte Carlo simulation. 
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A Determination of the D° Xe^v cocktail 

In this appendix, the exclusive semi-leptonic branching fractions, and a list of 
their averages are presented. This list, which is referred to as the Xe^u 
cocktail, is used to calculate the electron detection efficiency. The — > 
Xe+z/ cocktail is determined using world averages to obtain the following 
ratios: 

Rk* = BiD"^ ^ K*-e+iy)/B{D° ^ K-e+u) (16) 
R^ = B{D° 7r-e+z/)/^(D° ^ R-e+u) (17) 
Rp = B{D° ^ p-e+u)/B{D^ ^ K-e+u). (18) 

Experimental upper limits are used to obtain estimates for the unobserved 
modes: 

Rk(i270) = -B(D° ->Xf(1270)e+z/)/B(D° ^ir"e+z/) (19) 
i?K*(i430) = B{D^ ^ K*-{im)e+iy)/B{D'^ ^ R-e+iy). (20) 

The central value used for these unobserved modes is set to half the 90% 
confidence level upper limit and with an error equal to ±100% of the central 
value. 

The ratio of an exclusive channel to the inclusive rate is then obtained 
from the following formulas: 



S — 1 + Rk" + Rtt + Rp + Rk{1270) + Rk* (1430) (21) 

Xk = l/S (22) 

Xk* = Rr* I S (23) 

= RJS (24) 

Xp = Rp/S (25) 
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-^ii'(1270) = Rk{1270)/S (26) 
-^X*(1430) = Rk*{U30)/S. (27) 

Throughout this appendix the results are written in terms of the 
branching fractions. Results from the sector are converted into 
equivalent branching fractions using isospin and the measured and 
lifetimes. Also semi-muonic measurements are converted into semi-electronic 
results by correcting for the phase space difference between the muonic and 
electronic modes . In several of the tables, two averages are presented, one 
which includes all the data presented in the table, and another with CLEO 
results excluded. This is done to avoid double weighting in the CLEO data 
when performing calculations. 



A.l Rk* = B{D^ K*-e+iy)/ B{D^ R-e+u) 

There are two methods to measure this ratio: direct and indirect. The direct 
measurements, given in Table |^, can only be performed when both the K 
and K* modes are reconstructed through the same parent species within 
the same experiment. The indirect measurement compares the K*e^i' width 
measured in decays to the K~e^v width measured in decays, via 

Rk* = ^:r7^. : ^ X 



13{D+ K-TT+TT- 

B{D^ K-TX+) 
B{DO K^e+u) ^ 
B{D+ ^ ;^-7r+7r+) 
B{D^ ^ K-^+) Td+ ■ ^ ' 

Table contains the world average for B{D^ — K~ e^v) /B{D^ K^tt'^) 
and Table ^ contains the world average for B{D^ K*^e^v)/B{D^ 
K~TT'^7f^) where the CLEO measurements have been specifically excluded 
as these measurements are used in the direct determination of Rk*- To 
determine i?^?*"^^*, the ratio of normalizing modes Kmr/Kir presented in 
Table ITU] is used. Using the world average for this ratio of branching fractions 



and the / lifetime ratio p| the value for i?^?**'^^'^* is measured to be 
0.559 ± 0.068. Averaging Rfl^"^ and i?}^?^^^* yields 

Rj^. = 0.579 ± 0.049. (29) 
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Experiment Reference K*e^v/Ke^v 



CLE093 
CLE091 


i 

16] 


0.62 ±0.08 
0.51 ±0.19 


Average 






0.60 ±0.07 



Table 7: Direct measurements of the B{D K*e^v)/B{D Ke^v) ratio 
and their weighted average. 



Experiment 



Reference K e^v/K n 



CLE093 
E687 (94) 
CLE091 
E691 
E687 (90) 



113 

im 

\TE\ 

m 



0.978 ±0.052 
0.865 ±0.051 
0.86 ±0.07 
0.91 ±0.13 
0.84 ±0.19 



Average without CLEO 
Average 



0.869 ± 0.046 
0.906 ± 0.031 



Table 8: Measurements of the B{D^ K-e+u)/B{D^ K-ti+) ratio and 
their weighted average. The average without CLEO measurements is also cal- 
culated separately to avoid multiple use of the CLEO results in determining 
Rk*- 
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Experiment Reference K*^e^v/K tt+tt^ 



E691 

E687 

CLEO 

E653 

Argus 

WA82 



average 
average 
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4 
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m 
m 



without CLEO 



0.49 ±0.06 
0.59 ±0.07 
0.67 ±0.11 
0.48 ±0.11 
0.55 ±0.13 
0.62 ±0.17 



0.527 ±0.041 
0.547 ±0.038 



Table 9: Measurements of the B{D~^ 
and their weighted average. 



K*^e+u)/l3{D+ K- 



ratio 



Experiment Reference 



(%) 



K TT + TT 
(%) 



B(D+-^X-7r+7r+) 
B{D°^K-ir+) 



CLEO 
ARGUS 
ALEPH 
Mark III 
Mark II 
ARGUS 
HRS 
Mark I 
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26 
PI 



|28 



m 



3.91 ±0.19 
3.41 ±0.30 
3.89 ±0.33 

4.2 ±0.6 
4.1 ±0.6 
4.5 ±0.7 

4.50 ±0.94 

4.3 ± 1.0 



9.3 ±1.0 



9.1 ± 1.4 
9.1 ± 1.9 



8.6 ±2.0 



2.35 ±0.23 



average 


without CLEO 


3. 


84 ±0.18 


8.98 ±0.98 


2.34 ±0.28 


average 




3. 


87±0.13 


9.1 ±0.7 


2.35 ±0.18 



Table 10: Measurements of the hadronic normalizing modes, 
K~rr+, D+ K-tt+tt+ and their ratio. The CLEO result on B{D+ 
K^'K^'K^)/B{D^ K^-K^) is a direct measurement of this ratio, and is not 
obtained by dividing the individual CLEO results. 
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Experiment Reference Mode vr e^v/K e+z/ 



CLEO 1 
Mark III | 
CLEO 1 


34 
35 


1 


TT e 


V 


0.103 ±0.041 
0.115 ±0.051 
0.17 ±0.06 


average 












0.121 ±0.028 



Table 11: Measurements of the BiD^ vr e^v)/B{D^ K e^v) ratio and 
their weighted average. 



A.2 = B{D^ 7T-e+iy)/{D^ R-e+u) 

The Cabibbo suppressed decay — > tt^c^z/ has been observed at Mark III. 
CLEO has made measurements of both the rt^e'^v and the 

TT^e^u decay chains. There is factor of two due to isospin that is needed to 
convert the TT^e^u measurement to a n^e^u branching fraction. 



The results are presented in Table 11 



A.3 Rp = B{D^ p-e+iy) /B{D^ R-e+v) 

Fermilab experiment E653 has published an observation of four D'^ — > p^fi'^u 
events based on a kinematic separation of the Cabibbo suppressed p^fi~^u 
signal from the more copious K*^p^v mode |^. They measure B{D~ 



p^p+v)/B{D+ K*^p+v) = 0.044l|j:[Ji ± 0.014. To obtain Rp this mea- 
surement needs be corrected by the isospin factor and multiplied by Rk* 
which gives; Rp = B{D+ p^p+v)/B{D+ K*^p+v) x Rk* x Ip = 
(0.044t|]:[|i ± 0.014) X (0.579 ± 0.049) x 2 = 0.051 ± 0.037. For Monte Carlo 
generation it is assumed that the form factor ratios for D° — > p~e^v decay 
are identical to that of the well measured —>■ K*~e^v decay. 

A.4 B{D^ {R*7rye^iy) upper limits 

Searches for higher K^*^ resonances and possible non-resonant contributions 
to D semi-leptonic decay have been performed by the fixed target experi- 
ments [1^. Although no evidence for these decays has been demonstrated 



we include K-{1270)e+u and K*-{U30)e+u in the Monte 
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Mode 






Y 






1.0 


0.555 ±0.024 




K*-e+u 


0.579 ± 0.049 


0.321 ±0.021 






0.121 ±0.028 


0.067 ±0.015 


D"^ 




0.051 ±0.037 


0.028 ±0.020 




/sTf (1270)e+z/ 


0.03 ±0.03 


0.017 ±0.016 




ir*-(1430)e+z/ 


0.02 ±0.02 


0.011 ±0.011 


Sum 




1.801 ±0.077 





Table 12: The world average or estimate of the ratio of exclusive channels 
relative to the K~e^v decay mode, Rm = B{D^ me'^v)/B{D^ 

K~e^v). The third column, is the ratio of the exclusive rate to the sum of 
the exclusive rates, = B{D^ me+z/)/Sum. 



Carlo simulation. The decays are generated unpolarized and with the fol- 
lowing strengths and errors, Rk(i270) = B{D^ /Cf (1270)e+z/)/B(Z^° 
K-e+u) = 0.03 ± 0.03 and Rx-iuso) = B{D^ K*~{U30)e+u)/B{D^ 
K~e^i>) = 0.02 ± 0.02. It is assumed that any non-resonant contribution to 
the inclusive rate will have a similar electron momentum spectrum distribu- 
tion as these higher order modes. 

A. 5 Calculation of the Xe+u cocktail 

Table |12| summarizes the relative rates, Rm (relative to K^e^v) ob- 



tained in the previous sections. The sum of these rates is then used to 
determine the ratio of each exclusive rate to the sum of all the exclusive 



rates as per Eqs. 27. Table 13 contains the efficiencies for these exclusive 



modes to pass the selection criteria. 

A. 6 Comparison of the inclusive rate to the sum of the 
exclusive measurements. 

One of the most frequent comparisons in the literature ||l], |^, [l^ is the sum of 
the observed exclusive channels to the measured inclusive rate. The method 
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p{n) 


e{K-e+u) 


e{K*-e+u) 


e(7r e+z/) 


e(p e+z/) 


e(A7(1270)e+z/) 


e(i^*-(1430)e+z/) 


MeV/c 


{%) 


(%) 


(%) 


(%) 


(%) 


(%) 


225-250 


40.4 


34.4 


42.4 


38.1 


20.4 


10.9 


250-275 


42.8 


36.3 


45.4 


39.5 


22.9 


11.7 


275-300 


45.6 


38.6 


47.6 


42.6 


20.6 


12.0 


300-325 


46.2 


40.4 


49.2 


43.8 


23.4 


13.0 


325-350 


48.6 


41.0 


51.0 


46.6 


27.2 


12.0 


350-375 


50.8 


44.1 


54.6 


43.7 


30.7 


14.6 


375-400 


51.9 


46.1 


56.6 


48.2 


29.4 


19.5 


400-425 


53.9 


45.2 


57.7 


57.4 


21.4 


34.3 



Table 13: Efficiencies for the exclusive decay channels used in the Xe^u 
cocktail. 



of comparing the inclusive measurement to the sum of the ratio of exclusive 
measurements is presented here. 

The following set of equations are used to calculate the branching fraction 

for the observed exclusive decays: 

B{D^ K 



The sum of the observed exclusive rates is then: 

m 

{1 + Rk* + + Rp) (34) 

The quantities rff = B(D^ K-e+iy)/B{D^ K-7r+) and B(D° 
K~Tr'^) are common to all derived exclusive branching fractions, and thereby 
effect the entire scale. 





_ ^Ke+u 
— 'Kit 


X 


B{D'> - 


R- 


-7r+) 


(30) 




- r^e+u 

— ' Ktt 


X 


B{D^ - 


R- 


"7r+) X Rk* 


(31) 




— ' Ktt 


X 


B{D^- 


^R- 


"tt"*") X i?^ 


(32) 




— ^Ktt 


X 


B{D^ - 


^R- 


■7r+) X Rp 


(33) 



References 



29 



[1] Jeffery D. Richman and Patricia R. Burchat, UCSB-HEP-95-08 or 
Stanford-HEP-95-01, to be published in Reviews of Modern Physics. 
In this review of charm and beauty semi-leptonic decays, they obtain 
a value of [5.73 ± 0.25]% for the sum of the exclusive channels and 
[7.10±0.29]% for the inclusive D semi-leptonic decay branching fraction. 
This results in a difference between the inclusive rate and the exclusive 
rate of (19.3 ± 4.8)% of the inclusive rate. 

[2] L. Montanet et al. , Review of Particle Properties Phys. Rev. D50, 1994. 
Within the Note on Semi-Leptonic Decays of D and B Mesons, Part /by 
R.J. Morrison and J. D. Richman, they obtain [5.73±0.25]% for the sum 
of exclusive semi-leptonic D decay branching fractions and [7.01±0.62]% 
for the inclusive D semi-leptonic branching fraction. This results in a 
difference between the inclusive rate and the exclusive rate of (18.3 ± 
Vo of the inclusive rate. 



[3] D.S. Akerib et al. , (CLEO Collaboration), Phys. Rev. Lett. 71, 3070 
(1993). 

[4] A. Bean et al. , (CLEO Collaboration), Phys. Lett. B317, 647 (1993). 

[5] Y. Kubota et al. , (CLEO Collaboration), Nucl. Instr. Meth. A320, 66 
(1992). 

[6] Charge conjugation is implied throughout this paper, the word electron 
will be used to refer to both e~'s and e^'s. 

[7] Other sources of right and wrong sign electrons are — Z)° mixing and 
flavor changing neutral currents, — > e"'"e~. Experimental limits and 
theoretical expectations on these rare processes make their contribution 
to the right and wrong sign yields negligible 0. 

[8] N. Isgur, D. Scora, B. Grinstein and M. B. Wise, Phys. Rev. D39, 799 
(1989). 

[9] The D —> K*e'^v decay has three helicity states in the limit of negligible 



lepton mass. In Monte Carlo generation the world average pO|, |21|, ^ 
of the three form factors involved in the decay is used to determined the 
ratio of longitudinal to transverse alignments {Tl/'^t = 1-23 ± 0.13). 



30 



[10] When combining the CLEO result for B{D^ K-e+v)/B{D^ 
K-TT+) and the CLEO resuh for B{D° Xe+u)/B{D° K-tt+) we 
use the following values which have reduced systematic errors: 

B{D^ ^ K-e+u)/B{D'^ ^ K-7T+) = 0.978 ± 0.025 ± 0.032 (35) 
B{D^ ^ Xe+u)/B{D'^ ^ K-7r+) = 1.684 ± 0.049 ± 0.061. (36) 

[11] When combining the CLEO result for B{D^ K-e+u)/B{D^ 
K-TT+) and the CLEO resuh for B{D^ Xe+u)/B{D^ K-71+) we 
use the following values with reduced systematic errors: 

B{D^ ^ K"e+u)/B{D'^ ^ K^7r+) = 0.978 ± 0.025 ± 0.032 (37) 
B{D^ ^ Xe+u)/B{D'^ ^ K-7V+) = 1.684 ± 0.049 ± 0.061. (38) 

[12] This value for X^ includes the CLEO and non-CLEO results for 
R = B{D^ K-e+u)/B{D° K-tt+). The weighted average of 
R{CLEO) = 0.978 ± 0.041 (from [g) and R{non - CLEO) = 
0.869 ± 0.046 is R{wtd) = 0.930 ± 0.031. To calculate Xk one divides 
by the branching fraction ratio R' = B{D^ — > Xe^v)/{D^ K~7r^). 
Using only CLEO data to calculate Xk, the statistical and system- 
atic errors of R and R' are reduced because of the cancelation of the 
contributions from B{D^ K^n^), and one calculates Xk = R/R' = 
(0.978±0.041)/(1.684±0.078) = 0.581±0.036, where R and R' are taken 
from [1^. Using R{non — CLEO) to calculate Xk, one obtains Xk = 
R/R' = (0.869 ± 0.046)7(1.684 ± 0.109) = 0.516 ± 0.043. In this case 
the full errors for R' are used (see equation 9), because the non-CLEO 
measurement of R has no experimental correlation with the CLEO mea- 
surement of R'. To combine the CLEO and non-CLEO results for R to 
calculate Xk, one needs a weighted value of R' from the two values given 
above. R'{wtd) = 1.684 ± [0.078 f (correlated) + 0.109 f{uncor related)], 
where f{uncor related) = 1 - /{correlated) = (l/0.046)V[(l/0.046)2 + 
(1/0.041)2] = 0.438. This results in R'{wtd) = 1.684 ±0.092. Taking the 
ratio R{wtd)/R'{wtd) = (0.930 ± 0.030)/(1.684 ± 0.092) one calculates 
Xk = 0.552 ± 0.035. 

[13] K. Kodama et al. , (E653 Collaboration), Phys. Lett. B313, 260 (1993). 
P.L. Frabetti et al. , (E687 Collaboration), Phys. Lett. B307, 262 
(1993). 



31 



[14] Bacino, W. et al. , (DELCO collaboration) Phys. Rev. Lett. 43, 1073 

(1979). 

[15] M. Witherell in the Proceedings of the Lepton Photon Symposium, 
Ithaca, NY. 1993. 

John Cumalat in the Proceedings of the 7*'* Meeting of the American 
Physical Society Division of Particles and Fields, Fermilab 1992. 
D. M. Potter in the Proceedings of the Joint International Lepton- 
Photon Symposium and Europhysics Conference on High Energy 
Physics, Geneva, Switzerland, 1991 CMU-HEP91-19. 

[16] G. Crawford, et al. , (CLEO Collaboration), Phys. Rev. D44, 3394 

(1991) . 

[17] G. Bellini et al. , (E687 Collaboration), contributed paper to ICHEP, 
Glasgow, Scotland (1994). 

[18] J. C. Anjos et al. , (E691 Collaboration), Phys. Rev. Lett. 62, 1587 

(1989) . 

[19] P.L. Frabetti et al. , (E687 Collaboration), Phys. Lett. B315, 203 

(1992) . 

[20] J.C. Anjos et al. , (E691 Collaboration), Phys. Rev. Lett. 65, 2630 

(1990) . 

[21] P. L. Frabetti et al. , (E687 Collaboration), Phys. Lett. B307, 262 

(1993) . 

[22] K. Kodama et al. , (E653 Collaboration), Phys. Lett. B286, 187 (1992). 

[23] H. Albrecht et al. , (ARGUS Collaboration), Phys. Lett. B255, 634 

(1991) . 

[24] M. Adamovich et al. , (WA82 Collaboration), Phys. Lett. B268, 142 
(1991). 

[25] H. Albrecht et al. , (ARGUS Collaboration), Phys. Lett. B340, 125 

(1994) . 



32 



[26] The updated ALEPH result can be found in, P. Burchat, Review of 
Charm Physics, SCIPP-93-44. The original ALEPH result can be found 
in, D. Decamp et al. , Phys. Lett. B266, 218 (1991). 

[27] J. Adler et al. , (Mark III Collaboration), Phys. Rev. Lett. 60, 89 (1988). 

[28] R. H. Schindler et al. , (Mark II Collaboration), Phys. Rev. D24, 78 
(1981). 

[29] H. Albrecht et al. , (ARGUS Collaboration), Phys. Lett. B324, 249 
(1994). 

[30] S. Abachi et al. , (HRS Collaboration), Phys. Lett. B205, 411 (1988). 

[31] I. Peruzzi et al. , (Mark I Collaboration), Phys. Rev. Lett. 39, 1301 
(1977). 

[32] K. Kodama et al. , (E653 Collaboration), Phys. Lett. B286, 187 (1992). 

[33] F. Butler et al. , (CLEO Collaboration) CLNS 95-95-1324 Submitted to 
Phys. Rev. D. 

[34] J. Adler et al. , (Mark III Collaboration), Phys. Rev. Lett. 62, 1821 
(1989). 

[35] M. S. Alam et al. , (CLEO Collaboration), Phys. Rev. Lett. 71, 1311 
(1993). 

[36] K. Kodama et al. , (E653 Collaboration), Phys. Lett. B316, 455 (1993). 



33 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ex/95 11014vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ex/95 11014vl 



